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, Motion of test particles along rotating curved trajectories is considered. The problem is studied 

f"^ ■ both in the laboratory and the rotating frames of reference. It is assumed that the system rotates 

£SJ ' with the constant angular velocity uj — const. The solutions are found and analyzed for the case when 

» ! . the form of the trajectory is given by an Archimedes spiral. It is found that particles can reach infinity 

while they move along these trajectories and the physical interpretation of their behaviour is given. 
The analogy of this idealized study with the motion of particles along the curved rotating magnetic 
field lines in the pulsar magnetosphere is pointed out. We discuss further physical development (the 
t ' conserved total energy case, when ui ^ const) and astrophysical applications (the acceleration of 

, particles in active galactic nuclei) of this theory. 
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I. INTRODUCTION 



> 
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• Rotation and relativity are those two features of motion, which do not easily match with each other. Still in 
astrophysics, with its abundance of extremely strong electromagnetic and gravitational fields, there are situations 
where motion is both rotational and relativistic. Most prominent examples include swirling astrophysical jets in 
active galactic nuclei (AGNs) and quasars, innermost regions of black hole accretion disks, accretion columns in X-ray 
pulsars and plasma outflows in radio pulsar magnetospheres. In these kinematically complex astrophysical flows, where 
' rotation is interlaced with the relativistic motion of particles, the coexistence of these two features of the motion leads 
to observationally puzzling phenomena with sophisticated and ill-understood physical background. The interest to 
these flows is not new, but the upgrade of highly idealized models to more realistic, astrophysically relevant levels is 
| still related with major theoretical and computational difficulties. 

Some important and basic theoretical issues, related with the relativistic rotation, are not uniquely defined and 
often evoke controversial interpretations. One of the most notable examples is the "centrifugal force reversal" effect, 
originally found in [1], and later [2-4] studied in detail. It was argued that under certain conditions the centrifugal 
force attracts towards the rotation axis both for Schwarzschild and Kerr black holes. In Ernst spacetime, which 
represents the gravitational field of a mass embedded in a magnetic field, the centrifugal force acting on a particle 
in circular orbit was reported [5] to reverse its sign even twice! In the simplest case of the Schwarzschild spacetime, 
strictly and essentially speaking, it was found that below the radius of the spatially circular photon orbit an increase 
of the angular velocity of a test particle causes more attraction rather than additional centrifugal repulsion. This 
effect was interpreted by Abramowicz [4] in terms of the centrifugal force reversal - it was stated that in such cases 
the centrifugal force attracts towards the axis of rotation! 

This interpretation was criticized by de Felice [6,7] (see also [8]), who argued that the discovered effect could be 
attributed to the strength of the gravitational field and be explained in a way which preserves the repulsive character 
of the centrifugal force. The spirit of this approach — to save the intuitively appealing nature of the centrifugal force 
as of "something which pushes things away" [7] — is theoretically valid and practically convenient. After all, in general 
relativity, there is no implicit way to define the centrifugal force: in any case one needs to introduce some sort of 
"3+1" spacetime splitting and dub as the "centrifugal force" some Newtonian-like expression, which looks like it [7]. 
Moreover, de Felice found several interesting examples of the ambiguity of the global concept of 'outwards and pointed 
out at the deep interrelation of this problem with the definition of the centrifugal force in relativity. Abramowicz 
studied further the problem of the local and the global meaning of 'inwards and 'outwards [9] and showed that the 
centrifugal force always repels outwards in the local sense, while it may attract inwards, towards the centre of the 
circular motion, in the global sense! The theoretical scheme for the operationally unambiguous definition of the inward 
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direction was suggested by de Felice and Usseglio-Tomasset [10] and later this approach was used for the geometrical 
definition of the generalized centrifugal force [11]. 

Therefore the effect, discovered in [1], is indubitably a genuine relativistic effect, although its interpretation in terms 
of the "reversal" of the centrifugal force is not implicit and is largely the matter of definition. 

Same is true for another rotational effect, disclosed by Machabeli and Rogava [12], on the basis of the relatively 
simple and idealized special- relativistic gedanken experiment: motion of a bead within a rigidly and uniformly rotating 
massless linear pipe. It was shown that even if the starting velocity of the bead is nonrelativistic, after an initial 
phase of usual centrifugal acceleration, while the bead acquires high enough relativistic velocity, it starts to decelerate 
and after reaching the light cylinder changes the character of its motion from centrifugal to centripetal. It was found 
that when the initial velocity v > \/2/2 the motion of the bead is decelerative all the way from the pivot to the light 
cylinder. 

Certainly no real pipe may stay absolutely rigid, especially nearby the light cylinder. Besides, in order to maintain 
the uniform rotation of such a device, one needs an infinite amount of energy. Therefore the setup considered in [12] 
was highly idealized. The constant rotation rate assumption was replaced in [13] with a more realistic one: the total 
energy of the system "rotator+pipe+bead" was assumed to be constant. It was found that the moving bead acquires 
energy from the slowing down rotator, but under favorable conditions the bead deceleration still happens. 

The results of [12] were interpreted by its authors in terms of the centrifugal force reversal. De Felice disputed 
this interpretation [14] and argued that, also in this case, like in above-mentioned general-relativistic examples, the 
generalized definition of the centrifugal force may guarantee the absolutely repulsive character of the force. He pointed 
out that an inertial observer will never see the bead reaching the light cylinder, because all light signals from the bead 
are infinitely redshiftcd. It was also shown that the vanishing of the radial velocity of the bead at the light cylinder 
can be interpreted in terms of the corresponding vanishing of the bead's proper time. 

Despite the controversy of interpretations it is generally believed that rotational relativistic effects could operate in 
different astrophysical situations and might, hopefully, lead to detectable observational appearances. Recently, Heyl 
[15] suggested that the observed QPO frequency shifts in bursters are caused by a geometrical effect of the strong 
gravity, similar to the Abramowicz-Lasota centrifugal force "reversal" 1 . As regards Machabeli- Rogava gedanken 
experiment, it implies that radially constrained, relativistic and rotationally (centrifugally) driven motion shows 
inevitable radial deceleration near the light cylinder. Evidently this effect might occur in a number of astrophysical 
situations, where motion is constrained, rotational and strongly (special) relativistic. 

One of the most important class of astrophysical flows, where this effect could show up, is centrifugally driven 
outflows. In the context of pulsar emission theory they were first considered in the late 1960s by Gold [17,18] (for 
recent studies see e.g. Ref.19 and 20). For accreting black holes, both of galactic and extragalactic origin, Blandford 
and Payne [21] first noted that centrifugally driven outflows from accretion disks could be responsible for the launch 
of jets, if the poloidal field direction is inclined at an angle less than 60° to the radial direction 2 . 

Recently Gangadhara and Lesch [23] suggested that centrifugal acceleration, taking place as a consequence of the 
bead-on-the-wire motion similar to the Machabeli- Rogava [12] gedanken experiment, may account for the acceleration 
of particles to very high energies by the centrifugal forces while they move along rotating magnetic filed lines of the 
rotating AGN magnetosphere. They claimed that the highly nonthermal, X-ray and 7-ray emission in AGNs arises via 
the Comptonization (inverse-Compton scattering) of ultraviolet photons by centrifugally accelerated electrons. The 
same processes was critically re-examined by Rieger and Mannheim [24] and it was found that the rotational energy 
gain of charged particles is efficient but substantially limited not only by the Comptonization but also by the effects 
of the relativistic Coriolis forces. The specific nature of the propagation of electromagnetic radiation in the rotating 
frame of reference [14,25] is another aspect of this problem, which still needs to be taken into account. 

The whole philosophy of the 'pipe-bead' (ot 'rotator-pipe-bead') gedanken experiments was to mimic the common 
situation in relativistic and rotating astrophysical flows, where the plasma particles are doomed to move along the 
field lines of governing magnetic fields. While we consider relatively small length scales, the shape of the field lines 
can be approximated as being straight. However, on larger length scales the curvature of the field lines turns out to 
be important for the physics of the plasma streams, which are guided by them. The natural question arises: how 
the motion of the bead changes when the pipe is curved? In other words, how the dynamics of particles, prescribed 
to move along the fixed trajectories, change when the shape of their involuntary tracks of motion is not straight!? 
Obviously, this is not only a mere theoretical curiosity, but the issue which might have a tangible practical importance. 



1 However later [16] it was found that the Heyl's calculations contained the sign error and the real effect could hardly account 
for the observed frequency shifts in the type I X-ray bursts. 
2 For a rapidly rotating Kerr black hole the critical angle can be as large as 90° [22]. 
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In astrophysical situations the role of the "pipes" is played by the magnetic field lines, and the latter are always curved. 
Therefore, it is clear that the study of the motion of test particles along prescribed curved rotating trajectories is 
a necessary and important step for the ultimate building of a physically meaningful model of centrifugally driven 
relativistic particle dynamics for rotating magnetospheres of pulsars and AGNs. 

It is the purpose of this paper to address the above stated issue. In particular, in the next section, we develop spccial- 
relativistic theory of the motion of centrifugally driven particles on fixed nonstraight trajectories. The formalism is 
developed both for the laboratory frame (LF) and for the frame of reference rotating with the system (rotating frame, 
or RF). Equations of motion are derived and solved numerically. The detailed study is given only for the case when 
the angular velocity of the rotation is constant. However, we also outline the formalism for the astrophysically more 
realistic case of the conservative 'rotator-pipe-bead' system with perceptible exchange of energy between the bead 
and the rotator, leading to the variability of the angular velocity of the whole system. In the final section of the paper 
we discuss the results, consider the directions and aims of the future study, suggest and discuss those astrophysical 
situations, where the obtained results could be useful for the clarification of puzzling observational appearances of 
related astronomical objects. 



II. MAIN CONSIDERATION 



The ideal two-dimensional system, which we are going to consider, consists of three basic parts: the device of the 
mass M and the moment of inertia X, rotating with the angular velocity u(t), hereafter referred as the rotator, the 
massless but absolutely rigid pipe steadily attached to the rotator; and the small bead of the mass m and the radius 
equal to the internal cross-section radius of the pipe. The bead is put inside the pipe and can slide along the pipe 
without a friction. Evidently, instead of the pipe-bead dichotomy, since we are considering the two-dimensional layout, 
one may think about the 'wire-on-bead analogy, which is sometimes used [24]. 

Contrary to the [12], where a straight pipe case was studied, now we let the pipe to be an arbitrarily flat curve, 
mathematically defined by: 

<p = ip(r), (la) 

with [dip/dr = <p'(r)]: 

$ ee rip'(r). (lb) 
The dynamics of the system may be studied basing on two alternative assumptions: 

1. It makes the task simpler to suppose that the kinetic energy of the rotator Em is huge and always Em 3> E m ; 
i.e., despite the exchange of the energy with the moving bead, Em stays practically constant. Hence, the angular 
velocity of the whole "rotator-pipe-bead" system (henceforth referred as the RPB system) stays constant: 

lo = const. (2a) 



In this case the first part of the triple RPB system (the rotator) continuously supplies the bead with energy 
and helps to keep the angular velocity of rotation constant. Therefore, the problem reduces to the study of 
the pipe-bead double system (the PB system) with the constant rotation rate. In the case of the straight pipe 
(ip = ipo) the problem has exact analytic solution, found and analyzed in [12]. 

2. It is more realistic to assume that the rotator energy E M is finite, so that the whole RPB system is conservative 
E tot ee E M + E m = const. There is a perceptible energy exchange between the rotator and the bead: both the 
energy Em and the angular momentum Lm of the rotator are variable and, consequently, the angular velocity 
of the rotation can not stay constant: 

uj 7^ const. (2b) 

The problem with the straight pipe and variable rotation rate u(t) has no analytic solution. It was studied 
numerically in [13]. 

With either (2a) or (2b) assumptions the pipe is always assumed to be the passive part of the system. In order to 
mimic a magnetic field line it is assumed to be massless, having no share in the energy and/or momentum balance of 
the whole system. Still the role of the pipe — as the dynamic link between the rotator and the bead — is significant: 
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it provides the prescribed "guiding" of the bead motion in the rotating frame of reference and makes the trajectory 
of the bead known in advance. 

In this paper the dynamics of the gedanken system is studied in detail only under the first, easier, assumption of 
the constant angular velocity. The rout to the solution of the problem under the second assumption is also given, but 
its full study needs separate consideration and will be published elsewhere. 

There are two natural frames of reference, in which the dynamics of this system could be studied. The first, inertial 
one, is the laboratory frame (LF), where the observer measures the angular velocity of the rotator (and the pipe) to 
be u>(t), while the angular velocity of the bead is equal to: 

Q(t)=u;(t)+ip'(r)v(t), (3) 

and the dynamics of the moving bead is governed by the pipe reaction force acting on it. Note that v(t) = dr/dt is 
the radial velocity of the bead relative to the LF. 

The second frame, rigidly attached to the rotator and rotating with it (hereafter referred as the rotating frame, 
or the RF), is non- inertial, but quite convenient for the inspection of the motion of the bead along the curved pipe. 
This original approach implies embodying of the form of the pipe into the metric of the rotating frame. It was used 
in [12] for the straight pipe case and proved to be quite efficient for the case when u)(t) is assumed to be a constant. 
On the contrary, as we shall see later, the LF treatment appears to be handier when the second (2b) approximation 
(E to t = const and L tot — const, rather than u{t) = const) is chosen. That is why it is important to consider the 
problem both in the LF and in the RF. 



A. Uniformly rotating PB system 



First, let us consider the problem in the laboratory frame of reference (LF) and ascertain that it admits full 
(numerical) solution of the associated initial value problem. Second, let us consider the same problem in the rotating 
frame of reference (RF). We shall see that when the (2a) assumption of the constancy of the rotation rate is used the 
latter approach is mathematically easier and provides fuller information about the dynamics of the system. 



1. LF treatment 



The most straightforward way to approach the problem is to consider it in the laboratory frame of reference, in 
which the spacetime is Minkowskian: 

ds 2 = -dT 2 + dX 2 + dY 2 = -dT 2 + dr 2 + r 2 d(f> 2 . (4) 

We use geometrical units, in which G = c = 1. Note that azimuthal angle <f>, as measured in the LF, is related with 
the azimuthal angle ip, measured in the RF, via the obvious expression: <f> = ip + Lot. The pipe reaction force F is 
the dynamic factor constraining the bead to move along the pipe. It is easy to see (from the 4-velocity normalization 
gafiU a U 13 — —1) that the Lorentz factor of the moving bead is: 

1 (t) = [l-r 2 tf -v 2 ]- 1 ' 2 . (5) 
The angle between the radius- vector of a point of the pipe and the tangent to the same point is given by the relation: 

a = arctan<J>, (6) 
and the components of the reaction force, acting in the radial and azimuthal directions, are 

F r = -|F|sina = --^== f |F|, (7a) 



F* = |F|cosa=-=L=|F|, (76) 
V 1 + $ 



respectively. 
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Defining the physical components of the bead relativistic momentum [m(t) = mo^(t)]: 

P r = mv, 



(8a) 



P = mrn, (86) 
we can write the two components of the equation of motion in the following way: 

P r -rtP ct> = F r , (9a) 

P + ilP r = F+. (96) 
Combining these equations we can, first, derive the equation: 

P r + $P + Q($P r - P ) = 0. (10a) 

It is easy to calculate that: 

Q = ip'v + ip"v 2 , (106) 

rh = m7 2 [(0 + rtp"v)rflv + (v + r 2 ip'Vl)v\\ (10c) 

and using these relations together with (8) we can easily derive the explicit equation for the radial acceleration of the 
bead: 

ruoQ — j 2 rv(<p' + uv)(fl + rip"v) 
r= ' ' ' 

where 

A=[l- w V+$ 2 ] 1/2 . (12) 

The Eq. (11) being of the form f = G(r,r) admits full numerical solution, as the standard initial value problem, 
providing the initial position of the bead, ro, its initial velocity, vo, and the shape of the pipe, tp(r), are specified. 

Defining the spatial vector of the 2-velocity v = (v, rCl), we can calculate the absolute value of the reaction force 
|F| using the equation [26]: 

m = F-v, (13a) 

which, in our case, leads to: 

^fgr~ r m = ru>\F\. (136) 

It is also easy to verify that the following quantity: 

* = m(t) - LurP v = to 7(1 - t 2 luQ) = const(t), (14) 

is the constant in time. This allows to find the solutions of the problem as functions of the specific value of this 
constant. In the next subsection we will see what is the physical meaning of this parameter - it turns out to be 
proportional to the proper energy of the moving bead in the RF. 

One important class of a possible shape of the curved trajectory is Archimedes spiral, given by the formula: 

ip(r) = ar, a — const. (15) 

In this case, since <p" = 0, from (11) it is easy to see that f ~ fl, while (14) implies that |F| ~ fl as well. Therefore, in 
the case of the Archimedes spiral trajectory, we can predict that the asymptotic behavior of the functions fl(t), v(t), 
v(t), and |F(i)| will be similar. 
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2. RF treatment 



We see the LF treatment allows to solve the problem and to obtain the complete information about the dynamics 
of the bead motion along the fixed nonstraight (curved) trajectories. However it is quite instructive and much more 
convenient to consider the same problem in the frame of reference, rotating with the pipe-bead system (rotating 
frame - RF). In order to do this we, first, need to switch from (4) to the frame, rotating with the angular velocity u). 
Employing the transformation of variables: 

T = t, (16a) 
X = rcoscf) — rcos(ip + ut), (166) 
Y = rsiruf) = rsin(tp + ut), (16c) 

we arrive to the metric: 

ds 2 == -(1 - u 2 r 2 )dt 2 + 2uj 2 dtdip + r 2 dip 2 + dr 2 . (17) 

For the straight pipe (ip = ipo) case (17) reduces to the metric ds 2 = —(1 — u 2 r 2 )dt 2 + dr 2 , which was basic metric 
for the [12] study. Now, for a curved pipe, defined by the equation (1), (17) reduces to the following form: 

ds 2 = -(1 - u?r 2 )dt 2 + 2uar<$>dtdr + (1 + <5> 2 )dr 2 . (18) 

For the resulting metric tensor 

., ,, _ /-(1-wV), \ no . 



cur®, 1 + $ 

we can easily find out that 

A = [— detig^)} 1 / 2 = (1 - Lu 2 r 2 + $ 2 ) 1/2 , (20) 

and, apparently it is the same function A, defined previously by (13). 

For this relatively simple, but nondiagonal, two-dimensional spacctime we can develop the "1 + 1" formalism. Doing 
so we follow as a blueprint the well-known "3+1 formalism, widely used in the physics of black holes [27-29]. Namely, 
we introduce definitions of the lapse function: 



A 1 -u 2 r 2 +$ 2 

■ (21) 



and the one-dimensional vector j3 with its only component: 

9JL = 

9rr l + $ 2 V ' 

Within this formalism (18) can be presented in the following way: 

ds 2 = -a 2 dt 2 + g rr (dr + (¥ dt) 2 . (23) 

Note that for the metric tensor (19) t is the cyclic coordinate and, moreover, in the RF the motion of the bead 
inside the pipe is geodesic - there are no external forces acting on it. Hence the proper energy of the bead, E, must 
be a conserved quantity. Employing the definition of the four velocity U a =dx a / dr we can write: 

E = -U t = -U*\gtt + 9tM = const. (24) 
On the other hand, the basic four-velocity normalization condition g a {)U a U^ — —1 requires 

U* = [-g tt - 1g tr v - g rr v 2 ] ~ 1/2 , (25a) 
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this equation, written explicitly, has the following form: 

X..21-V2 



U* = [1 - w 2 r 2 - 2ujr<S>v - (1 + $ 2 )^ 2 



(256) 



Recalling the expression (3) for the angular velocity of the bead fl(t), measured in the LF, and the definition (5) of 
the Lorentz factor j(t) in the same frame of reference we can easily see that: 



u l = {i-r 2 n 2 -v 2 ]- 1/2 = i(t). 



(25c) 



It is important to note that the conserved proper energy of the bead, E, defined by (24) may be written in terms 
of the Cl(t) and j(t) functions simply as: 



E = 7(i)[l - r 2 (t)cofl(t)} = const. 



(26) 



Taking time derivative of this relation and rearranging the terms we will finally arrive to exactly the same Eq. (11) 
for the radial acceleration of the bead r as in the LF treatment. Note also that "J" = moE. 

But the convenience of the RF treatment goes much further. From (24) and (25a) we can derive the explicit 
quadratic equation for the velocity: 



{g 2 r + E 2 g rr )v 2 + 2g tr (g tt + E 2 )v + g tt (g tt + E 2 ) = 0, 



with the obvious solution: 



v = r 



Vgtt + E 2 

(g 2 r + E 2 g„) 



-gtrVgtt + e 2 ± ea 



(27) 



(28) 



The "1 + 1" formalism helps to write equivalents of the same equations in a more elegant form. Namely, if we define 
the radial velocity 



and corresponding Lorentz factor: 

then, instead of (25c), we will simply have: 

while from the (24) we obtain: 



V r =-(v + p r ), 
a 



7=(1-V 2 )- 1/2 , 



E = ^[a-0-V)]. 



Instead of (27) we will have [V 2 = g rr V r V r = V r V r , (3 2 = g rr (3 r (3 r = (3 r (3 r \: 



{0 1 + E 2 )V 2 - 2a(0-V) + (a 2 - E 2 ) = 0, 



with the solution: 



yr = 



P 2 +E 2 



a(3 r ± Ea 



lE 2 +[3 2 -a 2 



(29) 

(30) 
(31) 
(32) 
(33) 

(34) 



Note that the RF Lorentz factor defined by (30) and the LF Lorentz factor, specified by (5) do not equal each other 

7(*) ^ (35) 

which is the manifestation of the obvious fact that Lorentz factor is not an invariant physical quantity. One can see 
that for the j(t) the following quadratic equation holds: 



(a 2 - /3 2 )Y - 2aEj + (E 2 + /3 2 ) = 0, 



(36) 
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with the following solution: 



7(*) 



a 2 -f3 2 



aE ± \(3\^ f3 2 + E 2 — a 2 



(37) 



Therefore, above developed theory allows us to look for the solution of the initial value problem for the bead moving 
along an arbitrarily curved pipe. The thorough consideration of different particular cases is beyond the scope of the 
present paper. Instead, we will give representative solutions for one of the simplest kinds of a spiral - Archimedes 
spiral - defined by (15). 

The scheme for the complete inspection of the problem for any given initial value problem is the following: First, 
one specifies the initial location of the bead (ro) and its initial radial velocity (vo). The values of the <^'(ro) and $(ro) 
are fixed as soon as we specify the form of the pipe (p(r). The initial values for the f2(0) and 7(0) are given by (3) 
and (5): 



tt =lj + ip'(r )v Q , 



(38) 



70 



r 2 n 2 



-1/2 



while the value of the bead proper energy, according to (26), is given as: 

E = 70 [1 - rowfio]- 



(39) 



(40) 



Working with the Eq.(28), as the first order ordinary differential equation for the radial position r(t) of the bead at 
any moment of time, we can subsequently calculate all other physical variables. On the Fig.l the set of solutions is 
given for the case of the Archimedes spiral with a = —5, rotating with the angular velocity lo — 2 for the bead, which 
initially was situated right over the pivot of the rotator (r = 0) and had initial radial velocity v =0.1. 

These plots tell us that in the limit of large distance from the rotator the value of the radial velocity tends to the 
asymptotic value: 



lim inf v(t) 



= —to/a, 



(41) 



which, in this case, is equal to = 0.4. 



Radial distance 



Radial velocity 



Radial acceleration 




FIG. 1. Graphs for the radial distance r(t), velocity v(t), acceleration a(t), the Lorentz factor j(t), angular velocity fl(t) and 
the absolute value of the reaction force |F(i)| for the rotationally (lo = 2) driven bead, moving on the Archimedes spiral with 
a = —5. ro = 0, vo = 0.1. 



The angular velocity of the bead in the LF tends to zero, as well as the absolute value of the pipe reaction force, 
implying that at infinity the bead asymptotically reaches the limit of the force- free motion. This limit is understandable 
also analytically, because from (3) and (28) we can see that: 



lo EVa 2 



(42) 
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n-*£HZ. (43) 



From these expressions it is clear that this regime is accessible iff the condition \a\ > u) holds! Otherwise, the 
particle is not able to reach the infinity. 

Since the shape of the function r(t) is almost linear it is instructive to make plots for the functions v(r) for different 
values of the initial radial velocity vq, but with all other parameters of the initial value problem being the same. On 
the Fig. 2 we plotted these functions for eight different values of the initial radial velocity. We see that when vq = 
the movement of the particle is force-free (geodesic) and uniform during the whole course of the motion. Physically it 
means that for this particular value of the vo the shape of the pipe follows the geodesic trajectory of the bead, in the 
RF, for the metric (17) on the rotating 2D disk, so the bead moves freely, without interacting with the walls of the 
pipe. When u < «oo, the particle moves with positive acceleration and asymptotically reaches the force- free regime 
in the infinity. While, when vo > the character of the motion is decelerative, but the force-free limit is reached, 
again, when the bead heads to infinity. 

One more example of the latter behavior, similar to the case shown on the Fig. 1, but plotted for the initial velocity 
v = 0.5 > Woo = 0.4 is given on the Fig. 3. Here we see that, unlike the case given on the Fig.l, the acceleration of the 
bead is negative all the time and it reaches zero "from below" , taking less and less negative values. While the angular 
velocity of the bead relative to the LF Q(t) is also negative from the beginning but its absolute value decreases and 
reaches the zero as the particle tends to the infinity. 



Radial Velocity 




0.2 - // 
O.I — / 

O ^— ' ' ' 1 ' ' ' 1 1 

O 2 4 6 8 10 12 14 16 18 20 

r 

FIG. 2. Graphs for the radial velocity v(r), when the initial value of the v(r) is taken to be: 0.01, 0.1, 0.3, 0.4, 0.5 (force-free 
value), 0.7, 0.9, 0.99. lo = 0.1, a = -0.2. 
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Radial distance Radial velocity Radial acceleration 




FIG. 3. Graphs for the radial distance r(t), velocity v(t), acceleration a(t), the Lorentz factor j(t), angular velocity fl(t) and 
the absolute value of the reaction force |F(t)| for the rotationally (w = 2) driven bead, moving on the Archimedes spiral with 
a = —5. ro = 0, vo = 0.5. 



B. Conserved energy case 



When the bead accelerates, it continuously takes energy from the rotator. So, if one needs to keep the rotation rate 
constant, one needs to supply the system with energy from outside. This is, certainly, less realistic setup than the 
assumption that the "rotator-pipe-bead" system is conservative, viz. its total energy E tot is constant. In this case, 
however, the bead acceleration can not be permanent, because asymptotically it extracts all energy from the rotator 
and reaches the regime: u>(t) — > 0, Em — > and E m — > E tot . Clearly, in this situation, it is more convenient to study 
the dynamics in the laboratory frame of reference (LF), in which the rotator and the pipe are rotating rigidly with 
the time-dependent angular velocity uj{t). As regards the bead, since the shape of the pipe is curved, its angular 
velocity relative to the LF is given by (3). Since the pipe is considered to be massless and absolutely rigid, it does not 
contribute any amount of energy and/or angular momentum to the total energy E to t and angular momentum L to t of 
the system. The rotator for simplicity is assumed to be a sphere of the radius R and the mass M having the inertia 
moment 

1 = -MR 2 , (44) 
5 

the energy 

E M = |^ 2 (t), (45) 

and the angular momentum 

L M =lu}(t). (46) 
Note that (44-46) are nonrelativistic expressions. If initially, at t = to, 

wofl«l, (47) 

then it will remain nonrelativistic during the whole course of the motion, because it is assumed that the bead constantly 
extracts energy from the rotator, while the latter slows down so that the angular velocity cu(t) is a monotonically 
decreasing function of time. The (47) condition seems to be valid for the known fastest rotators in the Nature — 
pulsars. For the Crab pulsar, for instance, R = 1.2 x 10 6 cto, iv — 190.4Hz and consequently ujqR/c ~ 7.6 x 10~ 3 . 
Even for the fastest millisecond pulsars uj^R/c < 0.25. This justifies the usage of nonrelativistic (44-46) expressions 
in our analysis. 
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The remaining part of the threefold system — the bead — is assumed to be of the rest mass ijiq. Its angular velocity 
and radial velocity relative to the LF, at any given moment of time, are J7(£) and v(t) = r, respectively. Even when 
the initial radial velocity of the bead is nonrelativistic (v <C 1), it is still necessary to write relativistic expressions for 
its energy and angular momentum, because the bead gains energy, accelerates and sooner or later its motion becomes 
relativistic. Therefore, its energy and angular momentum must be written as: 

E m = m{t) = may(t), (48) 

L m = m -/(t)r 2 (t)n(t) = m(t)r 2 (t)n(t), (49) 

where the LF Lorentz factor j(t) is defined by (5). 

The system "rotator-pipe-bead" is conservative, there is no energy inflow from outside. In this sense it principally 
differs from the one considered in the previous section, where either the external energy source was necessary to keep 
the rotation rate u> constant or the rotator was assumed to possess an infinite amount of energy. Now, since the 
system is conservative, its dynamics are governed by the conservation laws of its total energy Etot = Em + E m and 
total angular momentum L to t = Lm + L m : 

^uj 2 (t)+m ol (t) = E tot , (50) 

Xw(t) + m {)1 {t)r 2 Q,{t) = Ltot. (51) 

And the solution of the problem reduces to the solution of these equations, linked with (9), for two unknown 
functions of time r(t) and u)(t) for an arbitrary initial value problem: initial location of the bead ro = r(0) and the 
initial value of the rotation rate w = w(0) 01 the whole system. The detailed study of this problem is beyond the 
scope of this paper and will be given in a separate publication. 

III. CONCLUSION 

The purpose of the present paper was to study the dynamics of relativistic rotating particles with prescribed, 
curved trajectories of motion in the rotating frame of reference. The work is a natural generalization of the gedanken 
"pipe-bead" experiment considered Machabeli and Rogava [12]. In that paper the authors considered the case of the 
straight rotating pipe and they found out that when the velocity of the bead, driven by the rotation of the whole 
device and sliding along/within the pipe, is high enough the character of the motion changes from the accelerated 
to the decelerated one. In particular, it was found that when the bead starts moving from the pivot (r = 0) of the 
rotating pipe with initial velocity vo > V2/2, the motion is decelerative from the very beginning. 

In this paper we consider the motion of rotationally driven particles along flat trajectories of arbitrarily curved 
shape. The practical motivation for this approach and its importance are related with the following two facts: 

1. The 'pipe-bead' (or the 'bead-on-the-wire') gedanken experiment is considered as a model for the study of dy- 
namics of centrifugally driven relativistic particles in rotating magnetospheres, in various classes of astrophysical 
objects, like pulsars [17-18,12,20] and AGNs [21-24]. The role of "pipes" is played by the magnetic field lines. 

2. The shape of magnetic field lines is always curved. It implies that for the large-scale, global dynamics of charged 
particles — driven by centrifugal forces and moving along curved field lines of rotating magnetospheres — it is 
important to know what qualitative changes occur when the form of the field lines is not linear but curved. 

In this paper we studied this problem, on the level of the idealized gedanken experiment, both in the laboratory 
(LF) and in the rotating (RF) frames of reference. For the simple example of the Archimedes spiral we found that 
the dynamics of such particles may involve both accelerative and decelerative modes of motion. 

One important difference from the linear pipe case [12] is that for the case of a curved pipe the motion of the 
bead is not any more radially bounded: there exist regimes of motion when the bead may reach infinity. This result 
has simple physical explanation. For the case of the linear pipe, rotating with the constant angular velocity u>q, the 
natural limit of the radial motion was given by the light cylinder radius, defined as Rl = Uq 1 ■ Now, in the case 
of the curved pipe, even when it rotates with the same constant rate, the bead can slide in the azimuthal direction, 
following the curvature of the pipe and having a variable angular velocity Ct(t). It means that now the role of the 
effective light cylinder is played by R>L(t) = fi(i) -1 , and, hence, all those radial distances become accessible, where 



11 



r(t) < Ri,(t). Therefore, if both r(t) and -Rl(£) are monotonically increasing functions, but the former stays always 
smaller than the latter (evidently it was the case in above considered examples for the Archimedes spiral) then the 
bead can reach infinity. 

Moreover, we found that there are special solutions, which are force-free during the whole course of the motion. 
These are simply geodesies in the two-dimensional rotating metric (17). In the LF the motion of the bead in this 
case is radial, because its angular velocity f2(£) stays zero all the time and, correspondingly, the light cylinder is at 
the infinity from the very beginning. The form of the trajectory in the RF in this case, ip(r) = — (vo/uo)r, is simply 
that trace, which a free bead could leave on the surface of the rotating disk during the course of its geodesic motion. 
Intuitively it is evident that if the pipe has this particular form the bead slides within it freely, without interaction 
with the walls of the pipe. 

We considered only one, simple, subclass of spiral trajectories as the representative example of the solutions, but 
the developed theory may be used for the study of the dynamics of particles moving along arbitrarily shaped flat 
trajectories. It means that this approach may find wide applications to different astrophysical situations where rotation 
impels plasma particles to move along curved magnetic field lines. 

We also gave basic equations and outlined the scheme for the solution of the more general version of the same 
problem, where the angular velocity of the rotating system is not assumed to be constant. Instead, it is assumed that 
the system rotator-pipe-bead is conservative and the rotator is allowed to exchange perceptible portions of energy 
with the bead. 

The formalism developed and the results found in this paper suggest some important directions of the further 
research, which could be physically and astrophysically relevant. 

First, it seems worthy to consider the case of the charged bead, which naturally would radiate while performing 
its nonuniform motion along the curved pipe. The radiative energy losses and the change of the angular momentum 
of the bead due to radiation could affect the dynamics of the bead and the whole system, considered, again, to be 
conservative. In the astrophysical context it could be interesting to see how the radiation of the bead would appear 
for the distant observer. 

Second, it is quite natural to try to extend the analysis for the 3-D fixed trajectories of motion, considering a family 
of axisymmctric trajectory lines and imitating the structure of the pulsar dipole magnetic field. This could bring us 
at least one step closer to the understanding of the radiative processes in pulsar magnetospheres. 

Third, sooner or later, we should address the fluid (plasma) problem and try to see how a continuous stream of 
fluid particles would behave, moving along rotating curved trajectories. This will comprise one more step closer to the 
reality of the pulsar environment and could help to relate with each other total energy losses of a pulsar, estimated 
through its slowing rate, and its radiation losses and energy taken away by centrifugally driven plasma, forming 
eventually the pulsar wind. 

These efforts would, hopefully, bring us to the construction of the unified theory of the pulsar magnetospherc, 
where the inertial aspects of the particle dynamics would be taken into due account. One could then try to test the 
theory with the existing empirical (observational) data about the energy deposited by pulsars into their winds and 
the energy they lose via their radio emission. This way we could have a clue as of how important inertial processes 
(often unfairly neglected) are in the dynamics of pulsars. 

Similar problems can be addressed also in the context of the centrifugal acceleration of particles in the jets in AGNs 
[21-24]. 

Yet another promising field of application is the accretion of plasma on strongly magnetized neutron stars, which 
is believed to lead to the appearance of X-ray pulsars. The final stage of the accretion is dominated by the dipole 
magnetic field of the accreting neutron star. It is normally assumed that the motion of accreted plasma particles is 
guided to the magnetic poles of the star by the rotating family of polar magnetic field lines. So, in a certain sense, 
this problem is of the same kinematic nature as the one related with radio pulsars, except that this time plasma 
moves from the region outside of the light cylinder towards the star. Certainly here, again, radiative effects and 
collective plasma effects are essential, so the simple 'one-particle treatment can give only very approximate picture of 
the involved physical processes. But taking into account the plasma fluid effects and the role of the radiation on the 
dynamics of infalling plasma streams, one could try to show how important the rotational (inertial) processes are for 
the dynamics of the flows infalling on strongly magnetized neutron stars and what is the influence of these processes 
on the observational appearance of related X-ray sources. 
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